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Abstract : A first insight into the phylogenetic analysis of the Orychophragmus violaceus complex was constructed from com- 


plete sequences of nuclear ribosomal DNA internal transcribed spacer regions (IT%5.8S) and chloroplast gene matK .The 


ITS results strongly supported the division of the O. wolaceus complex into two major groups: one is composed of O. hu- 


pehensis, O. taibaiensis, the other consists of O. violaceus and O. diffusus . The matK phylogenetic trees showed O. dif- 


fusus and O. violaceus were the basal members of the complex . From previous karyotype study and the present results, we 


considerated that the groups of O. hupehensis and O. taibaiensis should not be allocated in O. violaceus . And two possible 


ancestral cytotypes were also confirmed by our phylogenetic analysis . Moreover, the biogeographic analysis implied a pres- 


umed center of recent differentiation in the region from upper to mid Changkiang River to Qin Mountain . 


Key words: Biogeography; ITS; matK; Orychophragmus violaceus; Phylogenetic relationship 


Orychophragmus violaceus (Linn .) O .E .Schulz, 
a member of the Brassicaceae family, is closely related 
to Brassica (Wu, 2001; Luo et al., 1995) . This an- 
nual herb is native only to China where it is labeled as 


wild weed or cultivated as a garden plant or a potential 


edible-oil crop (Luo et al., 1995) . 

Since O. violaceus displays a high variability, its 
taxonomy has long been a controversial subject, with its 
variety identification and species status frequently re- 


vised . Thus Schulz (1903) recognized one species with 
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five varieties, while Tan et al. (1998) identified three 
additional species: O. taibaiensis Tan & Zhao, O. d- 
iffusus Tan & Xu, and O. hupehensis (Pampanini) Tan 
& Zhang . Among them, O. hupehensis was ever de- 
fined as a variety of O. violaceus in Flora of China 
(Zhou, 1987) . Later, Al-Shehbaz and Yang (2000) 
concluded that the variation of O. violaceus did not 
seem to follow any geographical distribution pattern . 
Hence they re-amalgamated these (including three ad- 
ditional species proposed by Tan et al.(1998)) into 
the original species O. violaceus as a highly variable 
species without any infraspecific taxa . Taking all of 
those evidences together, we can infer that the problem 
of reclassifying O. violaceus s variations and three addit- 
ional species cannot be effectively resolved by applying 
conventional morphological taxonomy . In the present pa- 
per, we adopted such scientific names as O. violaceus 
complex encompassing O. hupehensis, O. taibaiensis 
and O. diffuses, all of which were proposed by Tan et 
al. (1998), as well as the original species O. violaceus 
in light of the present findings . 
Orychophragmus hupehensis, O. taibaiensis and 
O. diffusus differ significantly from ©. [violaceus in 
both morphology (especially in leaf morphology) and 
distribution . O. violaceus and O. diffusus are biennial 
herbs mostly with ovate leaves . O. violaceus features 
narrowly ovate or oblongate leaves without lobe and 
erect stems growing upward, a few branches at the 
base, and some amplexicaul leaves at the top (Luo et 
al., 1995) . O. diffusus grows with diffuse caules and 
many branches at the base . All of its leaves are lyrata 
with serrate edge; the terminal leaf is cordate or reni- 
form and the leaf tip is blunt (Tan ef al., 1998) . 
However, O. taibaiensis and O. hupehensis are annu- 
als and characterized by the uppermost stem leaves hav- 
ing approximate triangle, being petiolate and not am- 
plexicaul . O. taibaiensis is different from O. hu- 
pehensis in that its uppermost stem leaves have slightly 
pinnate cleavage, cordate or ovoid triangles with apicu- 
late or obtuse ends, whereas for O. hupehensis its stem 
is underdeveloped, with all leaves largely pinnated 
cleavage, inverse-triangle and apex acute to acumi- 


nate . They are not overlapping but neighboring in their 


distribution . O. taibaiensis, native to Shaanxi, Meix- 
ian, Mt . Taibai of China, grows at an elevation of 
1150 - 1470 m, a height typical of cold climate of 
north-western China (Tan et al., 1998). O. hu- 
pehensis is, however, generally confined to Danjiankou 
and Shiyan of Hubei, and grows at an elevation of 450 
m to 1000 m where it is relatively warm . 

A recent phylogenetic analysis based on plastid 
and nuclear molecular markers suggested that O. vio- 
laceus was a sister clade to Brassica ( Warwick and 
Sauder, 2005) . However, this only provided one ac- 
cession of O. violaceus . As little is known of the evolu- 
tionary history of this highly variable taxa, molecular 
phylogenetic analysis is crucial in evaluating relations 
among the chromosomal races by using the mafK region 
of the chloroplast genome; and among different local 
populations by using the more variable internal tran- 
scribed spacer (ITS) sequences of nuclear ribosomal 
DNA (nrDNA) . The usefulness of both markers in 
species-level systematics and phylogeny of angiosperms 
has been widely demonstrated in recent studies (Diane 
et al., 2002; Winkworth et al., 2002; Li et al., 
2004; Selvi et al., 2006) . Based on our previous re- 
search, the goals of this study are to examine the cir- 
cumscription of O. violaceus, to investigate its major 
lineages, and to explain the evolutionary relationships 


within the O. violaceus complex . 


1 Materials and methods 
1.1 Plant materials 


Fifteen natural populations were sampled across the main 
distribution of the O. violaceus complex . These representative 
populations were from Hubei, Jiangxi, Shaanxi, Sichuan, Shan- 
dong, Jiangsu and Zhejiang Provinces and the cities of Beijing 
and Shanghai . As O. diffusus is rare and has a small population 
size, we only found one population in Shanghai . The geographi- 
cal origins of accessions are given in Table 1 and all abbreviations 
of the populations mentioned in the text follow as Table 1. 
Voucher specimens were deposited in the Herbarium of Sichuan 
University . Silica-gel dried samples of leaf tissue of each popula- 
tion were prepared for molecular analyses . 


1.2 DNA traction 


Genomic DNA was extracted following a modified 2x CTAB 
protocol (Doyle and Doyle, 1987) using samples of tissue cut 
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from leaves . The concentration of DNA in the various samples 
was determined by measuring the absorbance at 260 nm (Axo ) of 
a ten-fold dilution of each sample . The quality of all DNA prepa- 
rations was checked by agarose gel electrophoresis (0.7% we) 
in TAE buffer (1 mmo. EDTA, 40 mmoÑ®Ł Tris-acetate) con- 
taining 1 U@mL of ethidium bromide by comparison with a known 
mass standard . 
1.3 PCR amplification and DNA sequencing 

Analysis of the matK region was done for O. diffusus and for 
one accession of O. violaceus, O. hupehensis and two populations 
of O. taibaiensis with different ploidy level . The two primers 
matKF (5 -TAATACCTTATTTIGACTIGTATCGCACTAT-3' ) and 
matKR (5 -CCAAATCATTAAGATAAAGAATATCCAAA TACC- 
7 ) used in this study were designed on the basis of the constant 
regions of the published sequences of Acorus calamus (GenBank 
NC007407) , Brassica napus (GenBank X63088), Arabidopsis 
thaliana (GenBank AP000423) and Disporopsis undulata ( Gen- 
Bank AB029767) resulting from the alignment . Analysis of ITS 
nuclear rDNA was done for O. diffusus, the two populations of 
O. taibaiensis, the three of O. hupehensis and nine accessions of 
O. violaceus representing the distribution range (Table 1); the 
primers ITS1-18S (O Kane et al., 1997) and ITS4 (White et 
al., 1990) were used . 

For both matK and ITS, PCR amplifications were performed 
in a total volume of 50u 1 of reaction buffer, 1.5 mmo. MgC} , 
20 pmol of each primer, 200 Umof. of each dNTP, 1 U of Tag 
DNA polymerase (Fermentas) and 10 ng of template DNA . Re- 
actions were performed in a Peltier Thermal Cycler ( Bio-RAD 
DNAEngine) and programmed for an initial denaturation step (3 
min at 94°C) followed by 30 cycles of 1 min at 94°C, 1 min at 
55°C, and 2 min at 72°C . The last cycle was followed by a final 
incubation of 10 min at 72°C . The samples were then stored at 
4°C . Subsequently, 51 of each amplification mixture was ana- 
lyzed by agarose gel (1% v@W) electrophoresis in TAE buffer 
containing 1 U@mL ethidium bromide . The PCR reactions were 
purified from excess salts and primer using the Qiagen QIAquick 
PCR Purification Kit (Realtimes) . Automated DNA sequencing 
was performed directly from the purified PCR products using Big- 
Dye Terminator v . 2 chemistry and an ABI 3730 XL DNA Se- 
quencer ( Applied Biosystems) by Beijing Genomics Institute . 
The ITS region was not cloned as there was no evidence of the 
presence of two or more ITS sequences in the same individual . 
1.4 Sequene alignment and analysis 

DNA sequences and overlapping fragments were assembled 
and edited using SeqMan of DNAStar software package . The ITS 
sequences were checked for orthology to sequences of O. vio- 
laceus (GenBank AY722506), Vella spinosa (AY722498) and 
Brassica rapa (AF531563) . The two species were then used as 


outgroups for cladogram construction based on their position in 
the phylogeny of the tribe (Warwick and Sauder, 2005) . The 
sequence boundaries between the two spacers (ITS1 and ITS2) 
and the coding regions (5.8S) of nrDNA were determined by 
comparison with the published O. violaceus sequence ( Warwick 
and Sauder, 2005) . 

Multiple alignments were automatically performed using 
CLUSTALW in MEGALIGN of DNA Star (DNAstar Madison, 
WI), and then further examined and slightly modified manually . 
The matK sequences were aligned with the sequences of Sinapis 
alba ( AB354277) and Brassica rapa ssp. campestris (AB354276) ; 
the two were then used as outgroup representatives in the matK 
tree in light of the close relationship to Orychophragmus . Be- 
cause different analysis methods are sensitive to different biases in 
the data set, Baum ef al. (1994) suggested that analyzing data 
with multiple algorithms is desirable and that clades consistently 
supported in different analyses might be considered more robust 
than those supported strongly by one search method but contra- 
dicted by another . Phylogenetic analyses for each matrix were 
carried out by Bioctrl package ( http @@mnh. scu. edu. cr 
Bioctrl.htm) using maximum parsimony (MP) , maximum likeli- 
hood (ML), and Bayesian inference (BI) methods in PAUP’ 
version 4.0b10 ( Swofford, 2003), PHYML version 2.4.3 
(Guindon and Gascuel, 2003), and MrBayes version 3.0b4 
(Ronquist and Huelsenbeck, 2003), respectively . For MP anal- 
yses, heuristic searches were conducted with 1000 replicates of 
random addition, one tree held at each step during stepwise addi- 
tion, tree-bisection-reconnection (TBR) branch swapping, Mul- 
Trees in effect, and steepest descent off . Gaps were treated as 
missing data, characters were equally weighted, and their states 
were unordered . Intemal branch support was estimated with 1000 
bootstrap replicates (Felsenstein, 1985) . For maximum likeli- 
hood (ML) analyses, the model of sequence evolution fitting the 
data best was determined with Modeltest, version 3.06 (Posada 
and Crandall, 1998) . GIR + I + G model best fit the data 
(based on the alignment without matrix) . It was applied along 
with the Modeltest- estimated parameters in heuristic searches with 
one addition sequence replicate and TBR branch swapping . Nod- 
al robustness on the ML tree was estimated by the nonparametric 
bootstrap (1000 replicates for ITS and 100 for matK) . Bayesian 
analyses were accomplished in MrBayes version 3 .0b4 using the 
best-fit models upon Akaike information criterion (AIC; Akaike, 
1974) by using Modeltest 3.7 (Posada and Crandall, 1998) . In 
Bayesian analyses, trees were generated by running four simulta- 
neous Metropolis- coupled Monte Carlo Markov (MCMC) chains 
and sampling one tree every 1000 generations for 1 , 000, 000 
starting with a random tree . The posterior probability (PP) was 


used to estimate nodal robustness . 
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Table 1 


Taxa and code 


Orychophragmus violaceus HS 


0 0 O0 U D0 0 


Location 


Huashan Mountain 


310 


Taxa, origins, vouchers and accessions of the O. violaceus complex investigated and other species used in this study 


O. violaceus NJ Nanjing 
O. violaceus TS Mount Taishan 
O. violaceus JJ Jiujiang, JiangXi 
O. violaceus BJ Beijing 
O. violaceus XH HangZhou 
O. violaceus CD Chengdu, Shuangliu 
O. violaceus CS ChaoShan, Yuhang 
O. violaceus SNJ Shennongjia 
O. hupehensis SY Wudang Mountain 
O. hupehensis DJ Danjiangkou 
O. hupehensis WD Shiyan 
O. taibaiensis JK Jiaokou, Taibai Mountain 
O. taibaiensis HP Haoping, Taibai Mountain 
O. diffusus SH Shanghai 
O. violaceus 
Diplotaxis brachycarpa 
Vella spinosa 
Sinapis alba 
Brassica rapa subsp. campestris 
Brassica rapa 
Acorus calamus 
Brassica napus 
Arabidopsis thaliana 
Note:“ * ” indicate accession numbers of matK gene . 


2 Results 
2.1 Sequence characteristics 

Sequences were deposited in GenBank-EMBL- 
DDB and can be retrieved using the numbers in Table 
1 . The ITS region of nrDNA comprising both ITS se- 
quences (ITS1 and ITS2) and the 5.8S rDNA was am- 
plified by PCR from 15 populations of the O. violaceus 
complex (Table 1) . The aligned matK sequences were 
1357 bp in length with no insertions or deletions . As 
matK is highly conservative in Brassiceae (Marcus et 
al., 2001), the number of nucleotide substitutions was 
limited . Ingroup variation was low, with only 19 vari- 
able positions (1.4% ) . In the phylogenetic analysis, 
1282 positions were constant (94.47%), only 26 sites 
were parsimony- informative (1.92%) .In contrast, the 
ITS region was quite variable among taxa. The ITS 
aligned sequence data set was 553 bp in length . In- 
group variation was remarkable, with 240 positions 
(42.55% ) being variable . Adding the two outgroups for 
the purposes of the phylogenetic analysis, 425 characters 
were constant, 63 parsimony-uninformative and 64 


(11.57% ) parsimony-informative . 


Alt .(m) Source or vouchers GB No . IT9@matK 
600 L.R . Zhou 0701 EU306542 
90 L.R . Zhou 0702 EU306548 
1000 - 1500 L.R . Zhou 0703 EU306552 
200 L.R . Zhou 0704 EU306546 
44 L.R . Zhou 0705 EU306544 
70 L.R . Zhou 0706 EU306554 
500 L.R . Zhou 0708 EU306551 
300 L.R . Zhou 0709 EU306541 
900 L.R . Zhou 0707 EU306550@£U306555 * 
450 L.R . Zhou 0712 EU306543 
400 L.R . Zhou 0713 EU306545 
1000 L.R. Zhou 0711 EU306553@EU30655 8" 
1150 L.R . Zhou 0714 EU3065474EU306556 * 
1210-1250 L.R . Zhou 0715 EU306540@EU54318 1“ 
80 L.R . Zhou 0710 EU306549@EU306557 ` 
Warwick et al. (2005) AY72 2506 
Warwick and Sauder (2005 ) AY722445 
Warwick et al. (2005) AY722498 
Lu et al. (2007) AB354277 * 
Lu et al. (2007) AB354276 * 
Warwick et al. (2005) AF531563 
Zhang et al. (2005) DQ008868 “ 
Handa et al. (1992) X63088 * 
Sato et al. (1999) AP000423 * 


2.2 matK analysis 

The 50% majority-rule consensus tree of 83 equally 

parsimonious trees and the ML phylogenetic tree, generat- 
ed by successive weighting analysis of the mafK subma- 
trix, are shown in Fig .1 . The MP tree had a length of 82 
steps, a consistency index (CI) of 0.9759 and a retention 
index (RI) of 0.9394 . The overall topology of the MP tree 
was essentially identical to that produced by the maxi- 
mum-likelihood (ML) and Bayesian analysis (Fig.2) . 
Differences between Bayesian, MP and ML trees con- 
cerned mainly the bootstrap values at some nodes, which 
showed slight vanation from tree-to-tree . For instance, 
the clade comprising the O. violaceus - O. hupehensis 
- O. taibaiensis complex was supported by bootstrap 
values of 100% , 94% and 95% in the Bayesian, ML 
and MP trees respectively . However, in the ML tree, 
O. hupehensis and the two populations of O. taibaiensis 
were not separated as the other two trees showed . 

The phylogenetic analysis based on matK showed 
that the O. violaceus complex formed a strongly supported 
clade (100% bootstrap values in all trees) . Despite low 
genetic divergence for matK, rooting of the tree with S. - 
alba and B. rapa ssp. campestris, resulted in O. diffusus 
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Cryvchophragmus tabaiensis x= 


HP( 4} 





Orvchophragmus tabaiensis x= 12 
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Fig .1 The phylogentic trees of matK of different cytotypes of the O. violaceus complex 


Left tree: The 50% majority-rule consensus tree of 83 equally parsimonious trees generated by matK sequences . Right tree: The maximum like- 


lihood tree from matK sequences from different cytotypes of the O. violaceus complex . Bootstrap values are shown on branches when > 50% . 


Abbreviations of O. taibaiensis accessions follow Table 1 . Base chromosome numbers and ploidy are indicated for ingroup and outgroup taxa . 
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Fig.2 The Bayesian trees obtained from matK sequences of different cytotypes in the O. violaceus complex 


Bayesian bipartition posterior probabilities are shown on the branches . Abbreviations of O. taibaiensis follow Table 1 . Base chromosome numbers 


are indicated for ingroup and outgroup taxa “ *” indicated Bayesian bipartition posterior probabilities <0.5 . 


as the earliest divergent lineage, which formed a well 
supported clade (95% BS in MP tree) sister to O. vio- 
laceus . O. violaceus was, in turn, moderately supported 
as the second basal member of the complex sister to 
O. taibaiensis and O. hupehensis . Additionally, the se- 
quence data set indicated that O. hupehensis differed from 
O. taibaiensis in only two deletions (one of 2 bp and one 
of 1 bp) and one 1-bp insertion . 
2.3 ITS analysis 

The consensus MP phylogenetic tree (L= 197, CI 
=0.7766, RI = 0.8027) and the Bayesian tree de- 


rived from IT%%5.8S sequences was shown with boot- 
strap values in Fig .3 and the ML tree with distribution 
of other characters was shown in Fig .4 . The topologi- 
cal identity is found in three phylogenetic trees in a 
majority of clades . Differences between MP, Bayesian 
and ML trees were just the values of some clades and a 
little difference in the clade B2 . Compared to Bayesian 
and MP trees, the positions of two populations SH and 
BJ in the clade B2 is transposed by the ML analysis, 
whereas MP tree and Bayesian tree had additional dif- 


ferences in the clade B2 . For example, Bayesian anal- 
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ysis exhibits paratactic relations between the popula- 
tions of SH, XH, CS, TS in the second branch of 
clade B2, while in MP tree, these populations are hy- 
potactic relations in the same branch . This may be re- 
sult from the different algorithms used by different 
methods . 

All constructed ITS phylogenetic trees congruously 
suggested that the O. violaceus complex was divided 
into two major clades, which are strongly supported 
(MPBP&MLBP&BPP: SMA .00 in cladeA; MPBP® 
MLBP@BPP: 6630.00 inclade B) . Clade A was 
composed of two very well supported subclades: clade 
Al containing a tetraploid and a diploid accessions of 
O. taibaiensis ( (MPBP@MLBP@BPP: IPRA .00) 
and clade A2 containing three populations of O. hu- 
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pehensis ( (MPBP@MLBP@BPP: 8381.00) . 

Being the sister group to clade A, the other major 
clade B consists of two subclades: clade B1 and B2. 
Clade B1, the weakly supported subclade (MPBP& 
MLBP@BPP: 444%) .94) , is made up of the popula- 
tions of HS, JJ and SNJ . In this clade, populations of 
SNJ and JJ clustered together with a very high bootstrap 
value ( (MPBP@MLBP&BPP: 964A .00) , confirm- 
ing a very close relationship between them . The other 
more strongly corroborated lineage ( MPBP&#MLBP& 
BPP: 99901 .00), clade B2, was formed by the re- 
maining populations of O. violaceus and O. diffusus . 
However, O. diffusus, which is represent the cytotype 
of chromosome base number x = 10, was nested in 
clade B2 . 
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Fig .3 The Bayesian trees (left) and the maximum parsimonious phylogenetic tree in MP analysis (right) 


generated by ITS sequences of the O. violaceus complex 


Bootstrap values and Bayesian bipartition posterior probabilities are shown on or below the branches . 


Abbreviations of populations are followed Table 1 
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Fig. 4 ML tree inferred from the ITS data showing the distribution of selected character states (1-6) 


The results of ML bootstrap analysis are shown above or below the branches . These characters are: 1, the chromosome numbers; 2, the terminal 


leaves’ shape (A , ovoid triangles with apiculate or obtuse ends; A , inverse triangles with apiculate ends; vx , 


narrowly ovate or oblongate; 


©, cordate or reniform); 3, number of the lateral lobes of the uppermost leaves; 4, the uppermost leaves are amplexicaul ( + ) or not ( - ); 


5, chroatographic points of flavonoides from the leaves in Rf . =0.24- 0.25 (DY = dark yellow; BYG = bright yellowish green;“ - “ = no 


point) ; 6, chroatographic points of flavonoides from the leaves in Rf . =0.36- 0.38 (PB = pale blue; the other codes followed as above) . 


3 Discussion 
3.1 Phylogendic rdationships within the O. vio- 
laceus complex 

The branches found in the analysis above were 
consistent based on molecular data . These results 
showed that the O. violaceus complex could be divided 
into two sister groups, which were strongly supported 
by three phylogenetic methods (MP, ML and Bayes- 
ian) . One major group was composed of population 
O. hupehensis and O. taibaiensis (clade A), and the 
other of accession O. diffusus and O. violaceus (clade 
B) .However, the two major groups (A and B) showed 
a paratactic phylogenic relationships in the ITS trees, 
whereas the phylogenic analysis based on matK indi- 
cated that O. hupehensis and O. taibaiensis were diver- 
gent later than O. diffusus and O. violaceus . This may 


be due to a slower evolutionary rate of matK than that of 


ITS .In our phylogenic analysis, the matK trees could 
distinguish group A and group B more clearly, whereas 
the ITS trees just proved them as two separate clades . 
O. hupehensis and O. taibaiensis were sister 
groups separated from O. violaceus, but due to their 
high bootstrap support, each one became a true clade 
in the ITS phylogenetic trees . Both ITS and matK data 
displayed comparatively close relationships between the 
two cytotypes, suggesting that they might be the same 
ancestor yet subsequently divergent in different evolu- 
tionary patterns almost at the same time . Significantly, 
two cytotypes of O. taibaiensis, the diploid (JK) and 
tetraploid (HP) accessions, clustered together with pa- 
This 
suggested that the tetraploid was probably formed di- 


ratactic relations both in ITS and matK trees . 


rectly from the same ancestor of the diploid in one sin- 


gle evolutionary step, but not derived from the diploid 
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with the similar morphology . 

It could also be inferred from the ITS cladograms 
that the populations of O. violaceus in clade B differen- 
tiated into two subgroups with the same chromosomal 
number (Zhou ef al., in press), though with wide ra- 
diation yet little morphological variation . For clade B1, 
the SNJ and JJ populations, which are two isolated 
populations in the upper and middle reaches of the 
Changkiang River, were strongly supported a very close 
relationship between them . This suggested a relatively 
recent disjunction . Additionally, 
showed O. diffusus fell inside O. violaceus in clade 4, 


but obviously, it was the earliest clade separated from 


our ITS analysis 


other cytotypes in the matK cladogram with very high 
bootstrap value (95% in MP) .One possible cause of 
this obscure relationship may be the absence of near 
relatives of the rare cytotypes, which resulted in the 
very small number of informative sites . Another cause 
may also be the different resolving power of the two 
genes . However, it could be inferred from matK data 
that there was a sharp distinction between O. violaceus 
and other cytotypes . 

3.2 Taxonomic implication 
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Given the presence of clear karyotype and other 
different characters (Fig .3), this division of the O. - 
violaceus complex was, not surprisingly, quite obvious . 
The present reconstruction indicated that the indepen- 
dence of O. hupehensis and O. taibaiensis from O. vio- 
laceus had been karyologically paralleled by a differen- 
tiation in chromosome number, karyotype asymmetry 
and heterochromatin content . O. hupehensis was char- 
acterized by karyotype of’ 12m+2m”" + 6sm + 2sm”” 
and a higher content of constitutive heterochromatin 
(about 45.5% of the total karyotype length) (Zhou et 
al. in press) . However, O. violaceus displayed a dif- 
ferent karyotype off 16m + 2m" + 4sm + 2sm” or 14m 
+ 2m" + 6sm + 2sm™” . The cytotype O. diffusus, 
which was clustered with O. violaceus, also showed a 
diploid complement of 2n = 2x = 20 (Zhou ef al., in 
press ) different from populations of group A . 

Moreover, this remarkable divergence between the 
two major groups (A and B) was also supported by 
morphology and phytochemistry traits (Zhao, unpul- 
ished) . Members of clade A was characterized by tri- 
angular shape in general, but populations in clade B 


were mostly ovate leaves (Fig .5) . Furthermore, the 
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Fig .5 The comparison of the leaves morphology of 15 populations of the O. violaceus complex 


For codes of the populations, see Table 1 . The members of group A are in clade A of the phylogentic trees and the members of group B are in clade B 
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morphological variations among them were stable even if 
naturalized to the same field and environment (personal 
observation) . This indicated that the variation of phe- 
notype was due to genetic rather than environmental 
factors . The thin-layer chromatogram of flavonoides in 
the leaves from different populations of the O. violaceus 
complex (Zhao, unpulished) displayed significant dif- 
ference between the two major groups . Therefore, these 
features and the results of our study supported the exclu- 
sion of O.hupehensis and O. taibaiensis from O. vio- 
laceus, rather than the inclusion of the three into one . 

Based on previous taxonomic analysis, the exclu- 
sion of O. hupehensis, O. taibaiensis and O. diffusus 
from O. violaceus was advocated by Tan et al. 
(1998 ) . From our karyotypical and phylogenic re- 
search work, we proposed that the groups of O. hu- 
pehensis and O. taibaiensis should not be allocated in 
O. violaceus . Their populations are mainly character- 
ized with the uppermost stem leaves morphology of ap- 
proximate triangle, being petiolate and not amplexi- 
caul, which are different from O. violaceus . Because 
the name of O. hupehensis is early than O. taibaiensis, 
we recommend to recovery the species of O. hupehensis 
which is include O. taibaiensis . To confirm this, much 
more sampling and further analyses, such as micromor- 
phology, hybridization or other genetic and molecular 
studies, should be performed within the complex or the 
genus Orychophragmus . 
3.3 Evoution and biogeography of the O. vio- 
laceus complex 

The phylogenetic study of the species group, com- 
bined with analysis of the recorded geographic distribu- 
tion, allowed a tentative reconstruction of the major 
features of their biogeographic history (Nikulina et al., 
2007) . The molecular phylogenetic analysis based on 
cpDNA matK and nrDNA ITS gene revealed a signifi- 
cant population differentiation pattern consistent with 
the results of karyotype analysis of the O. violaceus 
complex . 

The matK analysis showed that the cytotypes of x 
= 10 and x= 12 were the basal clades of the O. vio- 
laceus complex tree . This result revealed the earlier di- 


vergence of x= 10 and x= 12 and the derived origin of 


O. hupehensis and O. taibaiensis, a strong proof of the 
hypothesis of two possible ancestral bases concluded 
from Karyotypical analysis, the lowest base x = 10 and 
the mostly conserved base x= 12 in the O. violaceus 
complex (Zhou etal., in press) . Although there were 
few informative matK sites, what underlined our data 
was the fact that this phylogenetic relationship was 
strongly supported . This proved, as Li et al. (2004) 
pointed out, that given the more conservative nature of 
matK, their molecular characteristics may well reflect 
the phylogenetic signals . 

The map of the O. violaceus complex distribution 
(Fig .6) showed that HS, JJ and SNJ accessions of 
O. violaceus, the members of clade B1, had sympatry 
with the group A in the upper and middle reaches of the 
Changkiang River to Mt . Qinling . In these areas, the 
O. violaceus complex exhibited higher diversity than in 
other area: two base number (x=11, 12) and a tetra- 
ploidy, as well as some populations from different major 
groups . Therefore we presumed that these areas may be 
the recent differentiation center for the O. violaceus 
complex . HS, JJ and SNJ accessions of O. violaceus, 
which maked up clade B1 and was located in this cen- 
ter, had lower chromosomal asymmetry than the popula- 
tions of clade B2 in other areas except for little morphol- 
ogy variation between them . Accordingly, we hypothe- 
sized that O. violaceus possibly radiated from the center 
into farther areas . This geographic pattern may be a 
consequence of undirectional or balanced selection . 

Additionally, notable genetic variations in O. vio- 
laceus were found via the phylogenetic analysis . The 
abundance of haplotypes detected by ITS analysis might 
result from changed environments or geographical seg- 
regation of this species . All the results indicated that 
radialized colonization probably occurred when new 
ecotypes were produced to adjust to new habitats . Eco- 
logical differentiation was perceived as a key factor for 
maintaining cytotypes in an allopatric condition 
(Thompson and Lumaret, 1992; Levin, 2002; Hus- 
band, 2004) . Taking into account the large genetic di- 
vergence between the cryptic species and their pa- 
raphyletic relationships, we inferred that the O. vio- 


laceus complex might have experienced rapid expansion 
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Fig .6 Location map of 15 populations in O. violaceus complex and distribution ranges of different groups 


Red dashed line: range of group A; black line: range of group B; green dashed line: range of subgroup B1 


and diversification or an even longer period of evolu- 
tionary independence (lack of gene flow between their 
locations) . Historical and ecological factors may cont- 
ribute to its allopatric speciation in China . 

This is the first insight into O. violaceus phyloge- 
netic analysis . We offer the preliminary hypotheses 
which are hopefully to be verified with additional sam- 
pling of more taxa and more characteristics . Further- 
more, specific divisions within the complex require ex- 
tensive efforts in redefinition in consideration of prob- 


lematic positions of these cryptic species . 
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